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SUMMARY 

The  onset  and  growth  of  delaminations  in  unnotched  [±30/±30/90/90] 

s 

graphite-epoxy  laminates  is  described  quantitatively.  These  laminates, 

designed  to  delaminate  at  the  edges  under  tensile  loads,  were  tested  and 

analyzed.  Delamination  growth  and  stiffness  loss  were  monitored  nondestruc- 

tively.  Laminate  stiffness  decreased  linearly  with  delamination  size.  The 

strain  energy  release  rate,  G,  associated  with  delamination  growth,  was 

calculated  from  two  analyses.  A  critical  G  for  delamination  onset  was 

determined,  and  then  was  used  to  predict  the  onset  of  delaminations  in 

[+1*5  /— 1+5  /0  /QO  ]  (n=l,2,3)  laminates.  A  delamination  resistance  curve 

n  n  n  n  s 

(R-curve)  was  developed  to  characterize  the  observed  stable  delamination 
growth  under  quasi-static  loading.  A  power  law  correlation  between  G  and 
delamination  growth  rates  in  fatigue  was  established. 


INTRODUCTION 


A  commonly  observed  failure  mode  in  laminated  composite  materials  is 
delamination  between  the  composite  layers.  Delaminations  may  develop  during 
manufacture  due  to  incomplete  curing  or  the  introduction  of  a  foreign 
particle;  they  may  result  from  impact  damage;  or  they  may  result  from  the 
interlaminar  stresses  that  develop  at  stress-free  edges  or  discontinuities. 
Furthermore,  delaminations  may  grow  -under  cyclic  loading.  Delamination 
growth  redistributes  the  stresses  in  the  plies  of  a  laminate,  and  may 
influence  residual  stiffness,  residual  strength,  and  fatigue  life.  Hence, 
a  fatigue  analysis  for  composite  materials  should  take  into  account  the 
presence  and  growth  of  delaminations. 

A  bibliography  of  experimental  and  analytical  work  on  delamination  is 
contained  in  reference  [l].  One  of  the  most  promising  techniques  for 
characterizing  delamination  growth  is  based  on  the  rate  of  strain  energy 
released,  G,  with  delamination  growth  [l].  Previous  work  [2]  has  shown 
that  the  cyclic  growth  rate  of  debonds  between  the  metal  and  composite 
components  of  reinforced  panels  could  be  correlated  with  G.  Measured 
critical  G  values  have  been  used  in  sophisticated  analyses  [3,^]  to  pre¬ 
dict  the  onset  of  edge  delaminations  in  unnotched  composite  laminates. 

In  the  present  study,  a  simple  technique  was  developed,  employing 
strain  energy  release  rates  to  characterize  the  onset  and  growth  of 
delaminations  in  a  composite  laminate.  First,  the  damage  that  developed 
in  unnotched  [±30/±30/90/90]s  graphite-epoxy  laminates  under  static 
tension  loading  and  tension-tension  fatigue  loading  was  determined.  Next, 
stress  distributions  generated  from  a  finite  element  analysis  were 
correlated  with  the  observed  damage.  Then,  during  quasi-static  test 
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loads,  delamination  growth  and  stiffness  loss  were  monitored  nondestruc- 
tively  to  relate  laminate  stiffness  and  delamination  size.  The  resulting 
test  data  and  analysis  were  used  to  derive  a  closed-form  equation  for  the 
strain  energy  release  rate,  G,  associated  with  delamination  growth.  Next, 
a  critical  value  of  G  for  delamination  onset  was  determined.  It  was  then 
used  to  predict  the  onset  of  delamination  in  [+^5ri/-'45r>/0  /90  ]  (n=l,2,3) 

laminates.  A  delamination  resistance  curve  (R-curve)  was  developed  to 
characterize  the  observed  stable  delamination  growth  during  quasi-static 
loading.  Finally,  a  power  law  correlation  between  G  and  delamination 
growth  rates  in  fatigue  was  established. 

SYMBOLS 

extensional  stiffness  matrix 

delaminated  area 

area  of  interface  containing  a  delamination 
strip  delamination  size 
incremental  strip  delamination  size 
delamination  growth  rate  in  fatigue 
coupling  stiffness  matrix 

half-width  of  laminate  cross  section 

empirically  determined  coefficients 
bending  stiffness  matrix 


[A] 

A 

A* 

a 

Aa 

da 

dN 

[B] 

b 

c,3 

[D] 
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axial  stiffness  of  a  partially  delaminated  laminate 


axial  laminate  stiffness  calculated  from  laminated  plate  theory 


axial  stiffness  of  a  laminate  completely  delaminated  along  one  or 
more  interfaces 


axial  stiffness  of  the  ith  sublaminate  formed  by  a  delamination 


initial  tangent  modulus  of  an  undamaged  laminate 


E11}E22,E33  lamina 


M  rate  of  stiffness  change  with  delamination  area 

dA 

rate  of  stiffness  change  with  fatigue  cycles 


G11’G12’G13  lamana  shear  moduli 

G  strain  energy  release  rate  associated  with  delamination  growth 


GI»GII»GIII  strain  energy  release  rate  components  due  to  opening, 
in-plane  shear,  and  out  of  plane  shear  fracture  modes 


GFEM  values  of  Gj,  Gi;r,  calculated  from  finite  element 

analysis 


Gc  critical  strain  energy  release  rate  for  delamination  onset 


delamination  resistance 


G  maximum  strain  energy  release  rate  in  constant  strain 

max 


amplitude  fatigue  test 


4 


h 

K 

£ 

N 


n 


t 


U,V,¥ 


U,V,¥ 

V 

da 

dA 

dig 

dA 


x,y9z 


x,y,z 


x. 
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£ 


£ 


O 


ply  thickness 

saturation  spacing  of  cracks  in  90°  plies 
gage  length  used  to  measure  axial  displacements 
number  of  fatigue  cycles 
number  of  plies  in  a  laminate 

laminate  thickness 

thickness  of  i"^  sublaminate  formed  by  delamination 
displacements  in  x,y,z  directions 
displacement  functions  in  x9y5z  directions 

material  volume 


rate  of  strain  energy  released  as  flaw  extends 
rate  of  work  done  by  applied  load  as  flaw  extends 

cartesian  coordinates 


nodal  forces  in  x,y,z  directions 

first  element  of  the  inverse  extensional  stiffness  matrix 
(Aij-1)  (i,j  =  1,2,3) 

nominal  axial  strain 


uniform  axial  strain  assumed  in  finite  element  analysis 
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nominal  axial  strain  at  onset  of  delamination 


£max  maximum  cyclic  strain  level  in  fatigue  * 

V12’V135V23  lami-na  Poisson’s  ratio  ♦ 

ax  axial  stress  in  a  ply 

Oc  remote  axial  stress  applied  at  onset  of  delamination 

Oz  interlaminar  normal  stress  between  plies 

SPECIMENS  AND  APPARATUS 

Unnotched  [±30/±30/90/90]  ,  T300-5208  graphite-epoxy  laminates  were 

s 

tested  in  tension.  This  laminate  was  designed  to  have  relatively  high 
tensile  interlaminar  normal  stresses  at  the  edges  resulting  in  the  formation 
of  a  delamination  [5].  Figure  l(a)  shows  the  delamination  that  developed 
along  the  edge. 

The  specimens  were  254  ram  (10  in.)  long  by  38  mm  (1.5  in.)  wide.  These 
eleven-ply  laminates  had  an  average  ply  thickness  of  0.l4  mm  (.0054  in.). 

Specimens  were  tested  in  a  closed-loop  hydraulic  testing  machine.  The 
specimen  length  between  the  grips  was  180  mm  (7*0  in.). 

Also  shown  in  figure  l(a)  is  a  pair  of  linear  variable  differential 
transducers  (LVDT's),  that  were  mounted  on  the  specimen  to  measure  displace¬ 
ments  over  a  102  mm  (4.0  in.)  gage  length.  To  prevent  slippage,  a  fast 
drying  glue  was  applied  to  the  central  5.6  mm  (.22  in.)  portion  of  the 

LVDT  mounts  where  they  touched  the  specimen.  During  Mstrain-controlledM  4 

loading,  these  LVDT’s  were  used  as  the  feedback  device  in  the  closed  loop. 
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Dye-penetrant  enhanced  radiography  was  used  to  monitor  delamination 
growth  through  the  specimen  width.  Diiodobutane  (DIB),  a  dye  penetrant 
opaque  to  x-rays,  was  injected  along  the  delaminated  edge.  The  film  was 
placed  immediately  behind  the  specimen.  While  still  mounted  in  the  test 
machine,  specimens  were  exposed  to  x-rays  generated  for  five  seconds  at  l8  kV 
from  a  portable  point-source  unit  positioned  386  mm  (l5«5  in.)  away  from  the 
specimen.  The  radiographs  showed  the  location  of  the  delamination  front 
(fig.  1(b)).  The  dark  outline  in  the  center  of  the  picture  is  the  "shadow" 
(x-ray  image)  of  the  LYDT  rods  used  to  measure  displacements. 

DAMAGE  DEVELOPMENT 

The  same  type  of  damage  developed  during  both  quasi-static  tension 
and  constant-amplitude,  tension-tension  fatigue.  First,  a  few  isolated 
cracks  formed  in  the  90°  plies.  These  were  followed  almost  immediately  by 
the  formation  of  small  delaminations  along  the  edge,  as  seen  in  figure  2(a). 
When  the  delaminations  formed,  the  number  of  90°-ply  cracks  increased 
significantly  along  the  delaminated  length  of  the  specimen.  Most  of  the 
90°-ply  cracks,  which  appear  as  horizontal  lines  on  the  radiographs  in  fig¬ 
ure  2,  extended  beyond  the  delamination  front  outlined  in  the  x-ray  photo¬ 
graph.  Many  of  the  ply-cracks  immediately  extended  halfway  across  the 
specimen  width.  As  loading  continued,  additional  delaminations  formed  and 
joined  with  original  delaminations.  Delaminations  grew  much  more  rapidly 
along  the  length  of  the  specimen  than  across  the  width  (fig.  2(b)). 
Eventually,  two  delaminations,  one  on  each  side  of  the  specimen,  extended 
along  the  entire  specimen  length  between  the  grips  (fig.  2(c)),  after  which 
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the  delaminations  continued  to  grow  across  the  width.  Loading  was  terminated 
when  the  delamination  front  reached  the  shadow  of  the  LVDT  rods  (fig.  2(d)). 

To  illustrate  the  location  of  damage  through  the  thickness,  a  few 
acetate  tape  replicas  of  a  delaminated  edge  were  made  [6],  Figure  3  shows 
photographs  of  two  replicas  and  a  portion  of  the  delaminated  edge.  As  shown 
in  figure  3(a),  ply  cracks  extended  through  the  thickness  of  all  three 
interior  90°  plies*  As  shown  in  figure  3(h),  delaminations  formed  and  grew 
in  -30°/90°  interfaces,  typically  shifting  from  one  interface,  through 
90°  ply  cracks,  to  its  symmetric  -30°/90°  counterpart.  However,  delamina¬ 
tions  did  not  shift  interfaces  at  every  90°  ply  crack  encountered.  As  shown 
in  a  photograph  of  the  delaminated  edge,  figure  3(c),  interface  shifting  did 
not  occur  in  a  regular  pattern. 

Besides  formation  of  delaminations  at  the  -30/90  interfaces  and  cracks 
in  the  90°  plies,  an  occasional  angular  crack  formed  in  the  innermost  -30°  ply. 
The  replica  of  the  edge  in  figure  ^(a)  shows  two  angular  cracks  in  the  inner¬ 
most  -30°  plies  originating  at  a  90°  ply  crack  tip  and  creating  delaminations 
in  the  +30/-30  interfaces.  As  shown  in  the  radiograph  in  figure  i|(b),  the 
delaminations  in  the  +30/-30  interfaces  were  small  and  triangular  in  shape. 
These  +30/-30  delaminations  often  temporarily  arrested  initial  -30/90  inter¬ 
face  delamination  growth  along  the  length  of  the  edge.  However,  -30/90  inter¬ 
face  delaminations  eventually  joined  up  and  grew,  whereas  isolated 
+30/-30  delaminations  usually  remained  small* 

STRESS  ANALYSIS 

Two  approximate  analyses  were,  used  to  obtain  quantitative  predictions 
of  the  onset  and  growth  of  delaminations.  The  first  was  a  quasi-three- 
dimensional  stress  analysis  that  yielded  stress  distributions  and  strain 
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energy  release  rates.  The  second  was  a  simple  rule  of  mixtures  analysis 
that  was  used  along  with  laminated  plate  theory  to  calculate  stiffness  loss 
and  strain  energy  release  rates. 

Because  delaminations  form  in  unnotched  laminates  as  a  result  of  the 
interlaminar  stresses  that  develop  at  the  edge,  a  quasi-three-dimensional 
finite  element  analysis  [7]  was  performed.  The  finite  element  analysis  was 
used  to  calculate  stress  distributions  in  the  [±30/±30/90/90]s  laminate 
for  a  unit  axial  nominal  strain  (e  =  l).  Some  details  of  the  analysis  are 
described  in  Appendix  A. 

Figure  5  shows  that  the  through-thickness  distribution  of  the  inter¬ 
laminar  normal  stress,  Oz ,  calculated  at  the  edge  is  compressive  in  the 
outer  30°  plies  but  reaches  a  relatively  high  tensile  value  at  the  -30/90 
interface  and  throughout  the  90°  plies.  Also  shown  in  figure  5  is  the 
approximate  az  distribution  through  the  thickness  calculated  from 
laminated  plate  theory  and  an  assumed  stress  distribution  across  the 
width  [8].  This  plot  also  shows  the  highest  tensile  stresses  to  be  at 

the  -30/90  interface  and  within  the  90°  plies. 

Figure  6  shows  a  distribution  of  0^  across  the  specimen  width,  near 
the  edge,  at  the  -30/90  interface,  as  well  as  a  distribution  of  the  axial 
stress,  0  ,  in  the  adjacent  90°  ply.  Both  o  and  O  have  high  tensile 
values  at  the  edge. 

These  stress  distributions  showed  reasonable  correlation  with  the 
observed  damage  that  developed.  Indeed,  examining  0Z  and  interlaminar 
shear  stress  distributions  are  helpful  in  identifying  likely  delamination 
sites.  However,  interlaminar  stress  distributions  calculated  from  finite 
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element  analyses  vere 


not  useful  for  modeling  damage  growth  quantitatively 
because  the  magnitude  of  calculated  peak  stresses  at  the  edge  varied  with 
mesh  size.  Furthermore,  linear  elastic  analysis  suggests  that  the  int  r 
laminar  stresses  at  ply  interfaces  can  become  singular  at  the  edge  [7].  This 
singular  behavior  would  preclude  the  use  of  a  failure  criterion  based  on 
maximum  interlaminar  stress  values.  Therefore,  an  alternate  approach,  based 
upon  strain  energy  release  rates,  was  adopted  to  quantitatively  describe  the 
onset  and  growth  of  delaminations . 

STIFFNESS  LOSS 

In  many  composite  laminates,  stiffness  loss  may  reflect  delamination 
growth.  Furthermore,  the  rate  of  stiffness  loss  with  delamination  growth  can 
be  directly  related  to  strain  energy  release  rates.  Therefore,  analysis  and 
experiments  were  performed  to  correlate  laminate  stiffness  and  delamination 

size. 

Rule  of  Mixtures  Analysis 

To  analyze  stiffness  loss  due  to  delamination,  a  simple  rule  of  mixtures 
analyses,  along  with  laminated  plate  theory,  was  used.  First,  the  stiffness 
(tangent  modulus)  of  a  balanced,  symmetric  composite  laminate  (fig.  7a)  was 
calculated  from  laminate  theory  [9,10]  as 


ELAM 


(1) 


where  XX1  is  the  first  element  of  the  inverse  extensional  stiffness 
matrix.  A, V1  (i.J  -  1,2,3),  and  t  is  the  laminate  thickness.  Next 


assuming  a  complete  delamination  in  one  or  more  interfaces,  and  using  the 
rule  of  mixtures  assumption  that  the  sublaminates  formed  undergo  the  same 
axial  strain  (but  no  longer  have  the  same  transverse  strains),  results  in 


E*  = 


n 


i=l 


t 


(2) 


where 

E*  =  stiffness  of  a  laminate  completely  delaminated  along  one  or 
more  interfaces 

=  the  laminate  stiffness  of  the  ith  sublaminate  formed  by  the 
delamination 

t^  =  the  thickness  of  the  i^^1  sublaminate 

Although  equation  (2)  represents  a  two-dimensional  formulation,  E*  will 
depend  upon  which  interfaces  delaminate.  This,  in  turn,  determines  the 
stiffness,  Eis  and  thickness,  tis  of  each  new  sublaminate.  Hence,  equa¬ 
tion  (2)  is  sensitive  to  the  through-thickness  location  of  the  delamination. 
For  the  [±30/±30/90/90]s  laminate,  assuming  a  delamination  in  both  -30/90 
interfaces  (fig.  7b)  equation  (2)  becomes 

,»  =  8E(±30)g  +  3E(9q)3 
11 

Finally  an  equation  for  the  stiffness,  E,  of  a  partially  delaminated 
specimen  was  developed  using  the  rule  of  mixtures.  Equal-sized  delaminated 
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strips  were  assumed  to  exist  at  both  edges  of  the  laminate  (fig.  7c).  Then, 
by  assuming  the  laminated  and  delaminated  portions  of  the  specimen  act  as 
independent  components  loaded  in  parallel,  the  rule  of  mixtures  yields 


E  ■  (E*  “  elam)  f  +  elam 


(4) 


A  more  general  form  of  equation  (k)  may  be  developed  by  assuming  that  the 
relationship  between  laminate  stiffness  loss  and  delamination  size  can  be 
represented  by 


E  ~  ELAM  =  A_  (5) 

E*  -  elam  a* 


where 


A  =  delaminated  area 


A*  =  total  interfacial  area 


Rearranging  equation  (5)  yields 


E  "  (E*  -  elam) 


a,  .  F 

A*  LAM 


(6) 


Equation  (h)  is  a  special  case  of  equation  (6)  where 


a  _  A_ 
h  A* 
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Experiments 


To  verify  the  linear  relationship  between  stiffness  and  delamination 
size  implied  by  equations  (4)  and  (6),  four  quasi-static  tensile  tests  were 
conducted.  The  specimens  were  loaded  in  a  strain  controlled  mode  until  a 
delamination  formed.  Then,  the  specimens  were  unloaded  to  ten  percent  of 
the  peak  nominal  strain,  DIB  was  placed  on  the  specimen  edges,  and  an  x-ray 
photograph  was  taken.  Next,  the  specimens  were  reloaded  in  increments 
of  250  pm/m  above  the  previous  maximum  strain  level.  This  procedure  was 
repeated  until  the  specimen  was  almost  totally  delaminated. 

During  each  loading,  output  signals  of  the  two  LVDT 1 s  were  averaged, 
and  load  deflection  curves  were  plotted  on  an  X-Y  plotter.  The  initial 
linear  portion  of  each  plot  was  used  to  calculate  laminate  stiffness 
corresponding  to  the  damage  recorded  in  the  previous  x-ray  photograph 
(fig.  8).  Delaminated  areas  recorded  on  the  photographs  within  the  102  mm 
(b  in.)  gage  length  were  measured  with  a  planimeter.  To  minimize  data  reduc¬ 
tion  error,  each  delamination  was  traced  three  times  and  measured  areas  were 
averaged.  Then,  a  strip  delamination  size,  a,  having  equal  area  over  the 
LVDT  gage  length  as  the  measured  delamination,  was  calculated  (fig.  9). 

Figure  10  shows  a  plot  of  normalized  stiffness,  E/Eq,  as  a  function  of 
normalized  delamination  size  a/b.  A  least-squares  regression  line  for  the 
data  indicated  that  E*  =  0.7^2Eo,  where  EQ  is  the  initial  tangent  modulus 
measured.  Hence,  a  total  delamination  in  the  [±30/±30/90/90]s  laminate 
would  result  in  a  25.8  percent  reduction  in  laminate  stiffness.  The  data 
agreed  with  the  linear  rule  of  mixtures  equation  (4),  normalized  by  ELAm, 
where  E*  was  calculated  from  equation  (3)  and  sublaminate  stiffnesses  were 
calculated  using  equation  (l)  (see  appendix  B)  with  material  properties 
from  reference  [ll]. 
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Finally,  equations  (2)  and  (4)  were  used  iteratively  to  calculate 
laminate  stiffness  for  specimens  having  -30/90  interface  delaminations  and 
concurrent,  although  small,  +30/-30  interface  delaminations  (see  Appendix  C). 
The  contribution  of  90°  ply  cracks  to  laminate  stiffness  loss  was  also 
considered  (see  Appendix  C).  However,  the  net  effect  of  both  secondary 
mechanisms  (+30/-30  delaminations  and  ply  cracks)  on  stiffness  loss  was 
negligible  for  the  [±30/±30/90/90]g  laminate. 

STRAIN  ENERGY  RELEASE  RATE 

For  an  elastic  body  containing  a  planar  flaw  of  area  A,  the  strain 
energy  release  rate,  G,  is  the  difference  between  the  rate  of  work  done, 
dW/cLA,  and  the  rate  at  which  elastic  strain  energy  is  stored,  dh/dA,  as  the 
flaw  area  increases  [12],  i.e., 

G  =  —  -  —  (7) 

dA  dA 

Assuming  that  a  nominal  strain,  £,  is  sufficient  to  extend  the  flaw,  the 
work  term  vanishes.  Then,  if  U  is  expressed  as  a  product  of  the  strain- 
energy-density  and  volume  of  the  body,  V,  substituting  Hooke's  law  into 
equation  (7)  yields 

G  =  -1/  £  M  (8) 

2  dA 

where  is  the  rate  of  stiffness  change  as  the  flaw  extends. 

dA 
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In  this  study,  the  body  was  a  tensile-loaded,  unnotched  composite 
laminate,  containing  edge  delaminations.  The  strain  energy  release  rate 
associated  with  the  growth  of  edge  delaminations  can  be  calculated  by 
assuming  two  strip  delaminations  (fig.  9)  where 

1/  =  2b  At 
A*  =  2b  5 

(9) 

A  =  2  l  a 
dA  =  25,  da 


Then,  substituting  equations  (9)  into  equation  (8)  and  differentiating 
equation  (h)  yields 


G  = 


(  E  -  E* 
l  LAM 


(10) 


Equation  (10)  may  also  be  derived  for  an  arbitrary- shaped  delamination  by 
substituting  equation  (6)  into  equation  (8),  differentiating,  and  noting 
that  V  =  A*t .  Hence,  as  indicated  in  equation  (10),  the  strain  energy 
release  rate  associated  with  delamination  growth  is  independent  of  the 
delamination  size.  The  magnitude  of  G  depends  only  on  the  laminate  layup 
and  location  of  the  delaminated  interface(s)  (which  determine  E^^-  and  E*), 
the  nominal  strain,  e,  and  the  laminate  thickness,  t. 

Furthermore,  the  strain  energy  release  rate  (eq.  (10))  may  have  con¬ 
tributions  from  any  of  the  three  components  Gj,  of  Gjjj ,  corresponding 

to  the  opening,  in-plane  shear,  and  out-of-plane  shear  fracture  modes.  In 
addition,  near  the  edge,  G  may  deviate  from  the  value  predicted  by 
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equation  (10),  which  was  developed  using  laminated  plate  theory  and  the  rule 
of  mixtures.  Therefore,  a  virtual  crack  extension  technique  was  used  with 
the  quasi-three-dimensional  finite  element  analysis  to  calculate  Gj, 
and  GIZI  as  a  function  of  delamination  size  (see  Appendix  A).  As  shown 
in  figure  11,  the  finite  element  analysis  indicated  that  the  total  G, 
represented  by  Gj  plus  Gjj  (Gjjj  was  negligible),  reached  the  value  pre¬ 
dicted  from  equation  (10)  once  the  delamination  had  grown  a  very  small  dis¬ 
tance  in  from  the  edge. 

Because  G  increased  rapidly  with  "a"  near  the  edge,  a  small  delamina¬ 
tion  that  formed  (for  whatever  reason)  at  the  edge  would  be  expected  to 
undergo  rapid  initial  growth.  This  behavior  was  observed  in  the  quasi-static 
tension  tests  used  to  generate  stiffness  data.  As  soon  as  a  delamination 
was  detected,  the  loading  was  stopped.  The  formation  and  growth  of  the 
delamination  to  some  finite  size  appeared  to  be  nearly  instantaneous.  There¬ 
fore,  G  calculated  from  equation  (10)  at  the  nominal  strain  where  delamina— 
tion  was  first  detected  was  considered  to  be  the  critical  value,  Gc ,  required 
to  form  the  delamination.  This  GQ  was  then  used  to  predict  the  onset  of 
delamination  in  other  laminates. 

DELAMINATION  ONSET 

To  predict  the  onset  of  delamination  in  other  laminates,  several  things 
were  done.  First,  tension  tests  of  [±30/±30/90/90]s  laminates  were  run  to 
determine  the  nominal  strain  level,  ec,  at  which  delamination  begins. 

Next,  ec  was  used  in  equation  (10)  to  predict  a  critical  Gc  for  the  onset 
of  delamination *  Then,  Gc  was  used  to  predict  the  nominal  strain  at  the 
onset  of  delamination  in  other  laminates.  A  more  detailed  description  of 
the  procedure  follows. 
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Critical  Gc  Determination 


First,  eighteen  [±30/±30/90/90 ]  graphite  epoxy  laminates  were  loaded 
monotonically  in  tension  at  a  rate  of  44.5  N/sec  (10  lbs/sec)  until  a 
delamination  was  detected.  The  load  level  corresponding  to  delamination  onset 
was  recorded  and  the  corresponding  applied  stress,  a  ,  was  calculated.  Then, 
to  determine  the  nominal  strain  at  the  onset  of  delamination,  e  ,  a  was 
divided  by  Ey^j,  calculated  from  laminated  plate  theory  using  the  following 
elastic  properties  from  reference  [ll] 

E11  =  138  GPa  (20-°  Msi) 

E22  =  15  GPa  (2.1  Msi) 

G12  =  5.9  GPa  (0.85  Msi) 
v12  -  °-21 

The  average  value  "was  3^70  ym/m*  In  addition,  was  determined  from 

LVDT  measurements  on  the  four  tests  conducted  to  generate  stiffness  data. 

In  each  of  these  four  tests,  the  load  deflection  plot  was  linear  until  the 
delamination  formed.  The  average  value  of  £c  where  the  load-deflection 
curve  deviated  from  linear  for  these  four  tests  was  also  3^+70  ym/m. 

Next,  £  was  substituted  into  equation  (10)  to  determine  G  •  Stiff- 
ness  E and  E*  were  calculated  from  equations  (l)  and  (3),  respectively, 
using  elastic  properties  from  reference  [ll].  The  average  laminate  thick¬ 
ness,  measured  with  micrometers,  was  1.51  mm  (0.059^  in.).  A  value 
of  137  J/m2  (O.78  in-lbs/in2)  was  calculated  for  the  critical  strain  energy 
release  rate. 
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Delamination  Onset  Prediction 


To  predict  the  onset  of  delamination  in  other  laminates,  equation  (10) 
was  inverted  to  yield 


2  G„ 


ec  = 


'^LAM  '  E*) 


(11) 


The  critical  G  determined  from  [±30/±30/90/90 ]  laminate  data  was  used 
c  s 

in  equation  (ll)  to  predict  the  nominal  strain  at  the  onset  of  delamination 
in  [+1%M5n/0n/90n]s  (n  =  1,2,3)  T300-5208  graphite -epoxy  laminates 
having  the  same  stacking  sequence  hut  different  thicknesses.  To  evaluate  E*, 
delaminations  were  modeled  along  both  0/90  interfaces  where  high  tensile  O z 
stresses  were  anticipated  [3]  on  the  basis  of  the  approximate  analysis  of 
reference  [8].  Delaminations  have  been  observed  in  the  0/90  interface  of 
a  [ ±4 5/0/90 laminate  in  reference  [12],  Then,  stiffnesses  Ej-^  and  E* 
were  calculated  from  equations  (l)  and  (2),  respectively,  using  elastic 
properties  from  reference  [ll]. 

A  ply  thickness  of  0.15  mm  (0.0057  in.)  was  determined  in  reference  [3] 
for  the  8-ply  (n  =  l),  l6-ply  (n  =  2),  and  2Uply  (n  =  3)  laminates.  In 
figure  (12),  predictions  of  £c  were  compared  to  eQ  values  calculated  by 
dividing  measured  values  from  reference  [3]  by  Because  the  data 

from  reference  [3]  represents  the  average  of  only  two  or  three  tests  for 
each  stacking  sequence,  only  preliminary  conclusions  can  be  drawn.  Never¬ 
theless,  the  good  comparison  indicates  that  Gc  may  be  independent  of  the 
ply  orientations  that  make  up  the  delaminating  interface.  For  example, 

Gc  may  be  the  same  for  delamination  onset  in  the  -30/90  interfaces  of 
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[±30/±30/90/90]s  laminates  and  for  delamination  onset  in  the  0/90  inter¬ 
faces  of  the  [+45  /-45ri/0n/90  ]  (n  =  1,2,3)  laminates.  In  fact,  equa- 
tion  (ll)  indicates  that  the  dependence  of  ec  on  stacking  sequence  and  the 
ply  orientations  that  make  up  the  delaminating  interface  is  accounted  for 
in  the  (Ej-^  -  E*)  term. 

In  addition 5  both  equation  (ll)  and  the  data  indicate  that  delaminations 
will  form  at  a  lower  nominal  strain  in  thicker  laminates  of  identical  stacking 
sequences.  According  to  reference  [k]  interlaminar  stress  distributions 
calculated  from  elastic  analysis  will  be  identical  for  all  /0  / 90  ] 

II  ii  IX  XI  o 

laminates  subjected  to  the  same  applied  stress.  Therefore,  the  authors  con- 
eluded  that  a  failure  criterion,  if  based  on  critical  interlaminar  stresses, 
would  not  predict  the  thickness  dependence  of  delamination  onset.  Hence, 
the  strain  energy  release  rate  appears  to  be  the  most  useful  parameter  for 
quantitatively  predicing  the  onset  of  delamination  after  the  delamination- 
prone  interface(s)  have  been  identified  from  a  stress  analysis. 

DELAMINATION  GROWTH 

Quasi -Static  Tension 

Edge  delamination  has  been  observed  to  be  a  stable  fracture  process 
in  laminates  subjected  to  tension  loading  [l,5,lU],  Hence,  the  applied 
load  must  be  increased  to  force  the  delamination  to  grow.  Tension  tests  on 
the  [±30/±30/90/90]g  laminates  confirmed  this  observation.  Unstable  growth 
of  the  delamination  through  the  width  did  not  occur  before  the  laminate 
failed  (fractured  into  two  pieces).  In  the  four  quasi-static  tests  con¬ 
ducted  to  generate  stiffness  data,  delaminations  did  not  grow  after  the  mean 
applied  load  exceeded  13,350  N  (3000  lbs).  The  four  laminates  eventually 
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failed  at  a  mean  load  of  20,000  N  (4500  lbs).  Therefore,  the 
[±30/±30/90/90]  laminate  is  well  suited  to  studying  stable  delamination 
growth. 

The  stable  growth  of  flaws  can  be  characterized  using  the  crack  growth 
resistance  curve  (R-curve)  concept  of  fracture  mechanics  [ 15 ] *  Therefore, 
a  delamination  resistance  curve  was  constructed.  Strain  energy  release 
rates,  G,  were  calculated  from  equation  (10).  Because  G  does  not  depend 
on  delamination  size,  it  appears  as  a  horizontal  line  in  figure  13.  The 
three  horizontal  lines  shown  are  G  values  calculated  for  a  single  specimen 
at  three  successive  nominal  strain  levels.  In  addition,  the  delamination 
resistance,  GR,  was  calculated  using  the  maximum  nominal  strain  in  equa¬ 
tion  (10).  However,  G^  was  plotted  as  a  point  corresponding  to  the  size 
(see  fig.  9)  of  the  delamination  created  by  the  nominal  strain.  As  shown 
in  figure  13,  the  curve  formed  by  all  such  points,  generated  during  the 
four  [±30/130/90/90]^  quasi-static  tension  tests,  constitutes  the  delamina¬ 
tion  growth  resistance  curve  (R-curve)  for  the  graphite  epoxy  specimens 
tested.  The  critical  Gc  used  in  the  previous  section  to  predict  the  onset 
of  delamination  represents  the  first  value  of  GR,  i.e.,  the  first  point  on 
the  R-curve.  Hence,  the  R-curve  characterizes  the  laminates  resistance  to 
delamination  growth  under  tensile  loading. 

If  the  R-curve  is  independent  of  the  ply  orientations  that  make  up 
the  delaminated  interface,  like  Gc  appears  to  be  (see  previous  section), 
then  equation  (10)  and  the  R-curve  can  be  used  to  predict  the  growth  of 
delaminations  under  quasi-static  tension  in  other  laminates.  As  the 
loading  is  increased,  G  can  be  calculated  from  equation  (10)  and 
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compared  with  the  R-curve  to  predict  delamination  size.  Hence,  delamination 
size  as  a  function  of  applied  load  could  be  predicted  for  other  laminates. 

Fatigue 

ExpejUrmnt6 

Constant  amplitude,  tension-tension,  strain-controlled  fatigue  tests 

of  [±30/±30/90/90]  graphite  epoxy  specimens  were  conducted  at  a  frequency 

s 

of  10  Hertz  and  a  strain  ratio  of  0.2.  Specimens  were  loaded  slowly  in 
tension  until  a  small  delamination  appeared.  Laminate  stiffness  was  measured 
during  this  initial  loading.  Then,  the  specimen  was  unloaded  and  delamina¬ 
tion  size  was  recorded  using  DIB-enhanced  x-ray  photography.  Next,  the 
specimen  was  reloaded  to  the  mean  strain  and  the  stiffness  of  the  delaminated 
specimen  was  recorded.  Then,  constant-strain-amplitude  cyclic  loading  was 
applied.  The  cyclic  loading  was  interrupted  at  specified  intervals  to 
measure  delamination  size  and  static  stiffness.  Hence,  a  photographic 
record  of  delamination  size  and  a  record  of  static  stiffness  as  a  function  of 
load  cycles  were  accumulated.  Cyclic  loading  was  terminated  when  the 
delamination  had  grown  across  most  of  the  specimen  width  and  the  front  was 
obscured  by  the  x-ray  image  of  the  LVDT  rods  (fig.  2(d)). 

V lamination  g/wtoth  Mold  mzoAvJimmtk 

Figure  1^  shows  a  typical  plot  of  delamination  size  as  a  function  of 
load  cycles.  A  strip  delamination  size,  a,  was  calculated  from  the  delami¬ 
nated  area  (fig.  9)  measured  from  x-ray  photographs  using  a  planimeter. 

Three  separate  planimeter  tracings  of  each  delaminated  area  were  performed  to 
minimize  data  reduction  error.  As  indicated  in  figure  l4,  once  the  delamina¬ 
tion  had  grown  over  the  entire  length  of  the  specimen  edge,  a  constant  growth 
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rate,  da/dN,  was  eventually  achieved.  Table  1  lists  the  growth  rates, 
da/dN,  determined  from  least  squares  linear  regression  analysis  of  the  data 
for  each  fatigue  test. 

Figure  15  shows  a  typical  plot  of  static  stiffness  loss  as  a  function  of 
load  cycles.  This  plot  also  became  fairly  linear  once  the  delamination  had 
grown  away  from  the  edge.  Stiffness  degradation  rates,  dE/dN,  were  calcu¬ 
lated  from  least  squares  linear  regression  analysis  of  the  data  over  the 
same  cyclic  range  used  to  fit  da/dN  data.  Because  stiffness  was  found 
to  be  linearly  related  to  delamination  size,  differentiating  both  sides  of 
equation  (1+)  with  respect  to  number  of  cycles  (N)  allowed  an  alternate 
determination  of  delamination  growth  rates,  da/dN,  from  measured  stiffness 
degradation  rates,  dE/dN,  i.e., 


da  _  /  b  \  dE 
dN  °  (e«  -  E jjJ  dN  • 


(12) 


Therefore,  delamination  growth  rates  for  the  edge  delamination  specimens 
could  he  estimated  without  measuring  the  delamination  size  directly  (Table  l). 
Hence,  the  [±30/±30/90/90’]  specimens  used  in  this  investigation  should 
he  useful  for  generating  baseline  delamination  growth  data. 

Data  coAAeJtation  wJtfh  anati^i^ 

Four  fatigue  tests  were  conducted  at  each  of  three  maximum  cyclic 
strain  levels.  The  maximum  cyclic  strain  levels,  £max ,  chosen  were  3000, 
3250,  and  3500  ym/m.  The  maximum  strain  energy  release  rate,  Gmax,  was 
calculated  from  equation  (10)  using  £  .  Therefore,  a  constant  £max  test 

H1UX 

was  also  a  constant  G„,„^  test.  The  mean  of  the  four  delamination  growth 

max 
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rates,  da/dN,  was  determined  for  each  of  the  three  G  levels.  A  power 
curve  of  the  form  =  c®max^  vas  fit  to  the  three  mean  values  of  da/dN 
and  the  three  values  of  G^^  using  a  least  squares  routine.  An  excellent 
correlation  was  achieved  for  both  delamination  growth  rate  measurement  tech¬ 
niques  . 

Figure  l6  shows  —  as  a  function  of  Gmax  for  growth  rates  estimated 
directly  from  measured  delaminated  areas .  The  growth  rates  for  all  four 
tests  at  each  Gmax  level  are  shown  along  with  the  mean  growth  rate.  Fig¬ 
ure  17  compares  the  least-squares  power  law  fits  using  direct  area  measure¬ 
ments  and  indirect  stiffness  estimates  of  delamination  growth  rates.  Because 
both  da/dN  and  G  can  be  calculated  without  direct  measurement  of  delamina¬ 
tion  size,  the  stiffness  technique  presents  a  relatively  simple  means  of 
generating  data  to  determine  the  dependence  of  empirical  parameters  c  and  3 
on  differences  in  load  history,  frequency,  temperature,  etc. 

With  the  parameters  c  and  3  determined,  the  power  law  might  be 
applied  to  other  laminates.  The  calculated  strain  energy  release  rates 
could  be  used  with  c  and  3  measured  from  edge-delamination  baseline 
tests  to  determine  delamination  growth  rates. 

CONCLUDING  REMARKS 

Methodology  Summary 

A  methodology  for  analyzing  the  onset  and  growth  of  delaminations  in 
composite  materials  was  formulated  based  on  findings  in  the  current  work 
and  in  the  literature.  The  methodology  is  as  follows. 

First,  a  stress  analysis  of  the  particular  material,  configuration, 
and  loading  must  be  performed.  The  analysis  should  establish  where 
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delaminations  will  be  located.  For  unnotched  laminates ,  the  approximate  analyses 
discussed  in  this  paper  [8,1 6]  might  be  adequate.  After  the  location  is  established 
the  delamination  onset  and  growth  can  be  characterized  quantitatively  using 
strain  energy  release  rates.  Such  a  characterization  incorporates  the  influ¬ 
ence  of  material  volume.  Furthermore,  determination  of  the  singular  stress 
field  at  the  delamination  front  is  not  required. 

The  current  work  indicates  that  a  total  strain  energy  release  rate 
may  be  sufficient  to  characterize  the  onset  and  growth  of  edge  delaminations 
in  tensile-loaded  coupons .  More  work  is  required  to  determine  if  this  is 
true  for  other  configurations  and  loadings.  However,  if  strain  energy 
release  rates  must  be  separated  into  Gj,  G-q,  and  G.^  components  using 
a  numerical  analysis,  then  the  dependence  of  these  calculated  components 
on  grid  size  should  be  carefully  checked  and  documented. 

Immediate  Results 

A  simple  rule  of  mixtures  analysis,  using  laminated  plate  theory, 
indicated  that  laminate  stiffness  was  a  linear  function  of  delamination 
size.  The  analysis  accurately  predicted  stiffness  loss  due  to  edge 
delaminations  in  [±30/±30/90/90]s  graphite  epoxy  laminates.  The  linear 
stiffness  relationship  was  used  to  derive  a  closed-form  equation  for  the 
strain  energy  release  rate,  G,  associated  with  delamination  growth  in 
unnotched  laminates.  The  simple  G  equation  was  used  to  predict  delamina¬ 
tion  onset  in  [+fc5nM5n/0n/90n]s  (n  =  1,2,3)  laminates  using  a 

critical  G  determined  from  [±30/±30/90/90]s  laminates.  Stable 
c 

delamination  growth  in  the  [±30/±30/90/90]s  laminates  was  characterized 
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by  developing  a  delamination  resistance  curve  (R-curve).  Delamination 
growth  in  [ ±30/ ±30/ 90/90]  laminates  in  fatigue  was  characterized  by 

developing  a  power  law  correlation  between  G  and  delamination  growth  rates. 

Potential  Applications 

Preliminary  predictions  of  delamination  onset  in  [-«-1+5n/-i+5n/On/ 90n] g 

laminates  using  a  critical  Gc  determined  from  tests  on  [±30/±30/90/90]s 

laminates  indicated  that  Gc  may  be  independent  of  the  ply  orientations 

that  make  up  the  delaminating  interface.  If  this  is  true  for  stable 

delamination  growth  with  increased  tensile  load  and  delamination  growth 

in  fatigue,  then  the  delamination  resistance  curve  (R-curve)  and  power 

law  developed  on  [±30/±30/90/90]  laminates  can  be  used  to  predict 

s 

delamination  growth  in  other  laminates. 
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FINITE  ELEMENT  ANALYSIS 

Formulation 

The  quasi-three-dimensional  finite  element  analysis  was  developed 
in  reference  [7].  A  displacement  field  of  the  form 

u  =  £Qx  +  U(y,z) 

v  =  V(y ,z )  (Al) 

w  =  W(y ,z) 

was  assumed,  where  was  a  prescribed  uniform  axial  strain.  Eight-noded 

quadrilateral,  isoparametric  elements  with  three  degrees  of  freedom  per  node 
were  used  to  model  a  cross  section  along  the  specimen  length  (fig.  18).  Only 
one  quarter  of  the  cross  section  was  modeled  due  to  symmetry  conditions.  Each 
ply  was  modeled  with  one  element  through  its  thickness  except  for  the  cen¬ 
tral  90°  ply,  which  was  modeled  with  one  element  through  its  half-thickness. 
The  graphite  epoxy  unidirectional  properties  [ll]  used  in  the  analysis  were 

E11  =  138  GPa  (2°'°  Msi) 

E22  =  E33  =  15  GPa  ( 2 -1  Msi) 

( A2 ) 

g12  =  G13  =  G23  -  5.9  GPa  (0.85  Msi) 
v12  ~  v13  =  v23  =  °*21 
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Virtual  Crack  Extension  Technique 

Previously,  a  virtual  crack  extension  technique  has  been  applied  to 
finite  element  analysis  of  delaminations  [1,4].  In  this  technique,  the  work 
required  to  close  the  delamination,  expressed  in  terms  of  nodal  forces  and 
displacements,  is  assumed  to  be  equivalent  to  the  strain  energy  released  as 
the  delamination  extends  due  to  a  constant  nominal  strain.  This  technique 
greatly  simplifies  the  computation  of  Gj,  Gj  j  5  and  Gjjj  because 

knowledge  of  the  singular  stress  field  near  the  crack  tip  is  not  required. 

Figure  19  illustrates  the  technique  as  it  was  used  with  the  quasi-three- 
dimensional  finite  element  analysis.  First,  the  nodal  forces  were  calculated 
for  an  initial  delamination  of  size  a.  Then,  the  delamination  was  extended 
an  amount  Aa  and  the  resulting  nodal  displacements  at  the  same  location 
were  calculated.  The  expressions  for  Gj,  G-^,  and  Gjjj  in  terms  of 

nodal  forces  and  displacements  were 


where,  for  example,  Z1  represents  the  force  in  the  z  direction  at  node  b 

b 

calculated  from  element  1. 

Figure  11  shows  G^  and  Gjj 
Gm  was  negligible  for  this  case. 


as  functions  of  delamination  size. 

Strain  energy  release  rates,  G_„, 

FEM5 


27 


APPENDIX  A 


calculated  using  equations  (A3),  where  normalized  by  the  constant  value  pre¬ 
dicted  from  equation  (10).  The  total  G,  represented  by  Gj  plus  G^. j , 
reached  the  value ^  predicted  by  equation  (10)  once  the  delamination  had  grown 
a  very  small  distance  in  from  the  edge.  Furthermore,  the  total  G  calculated 
was  not  sensitive  to  mesh  refinement.  However,  the  values  of  Gj  and  Gjj 
calculated  were  sensitive  to  mesh  refinement. 

Figure  11  shows  the  nodal  discretization  used  in  the  finite  element 
analysis.  Four  different  initial  delamination  sizes,  a,  were  modeled.  For 
each  a,  the  delamination  was  grown  in  ten  increments,  Aa,  equal  to  one  tenth 

^Because  the  finite  element  analysis  models  only  one  quarter  of  the 
laminate  cross  section  due  to  symmetry,  the  analysis  assumes  four  delaminations 
grow  simultaneously,  two  on  each  side,  located  in  the  -30/90  interfaces.  How¬ 
ever,  in  the  derivation  of  equation  (10),  only  a  single  delamination  was  con¬ 
sidered  on  either  edge  to  be  consistant  with  the  physically  observed  behavior 
shown  in  figure  3(c)  and  illustrated  in  figure  20(a).  Making  the  same  assump¬ 
tions  as  the  finite  element  analysis  would  have  yielded 

V  =  2b£t 
A*  =  4b& 

A  =  kSLa, 
dA  =  da 

which  would  result  in  G  =  £-£-  -  E*J.  Hence,  the  values  calculated  from 

the  finite  element  analysis  using  equations  (A3)  were  doubled  to  be  consistent 
with  equation  (lO)  and  the  observed  physical  behavior. 
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of  the  initial  delamination  size.  As  a  was  increased,  the  mesh  refinement 
changed  accordingly.  For  the  smallest  delamination  size,  G^  continually 
increased  as  the  delamination  grew,  but  G^  was  constant  for  a.l  i  ten  incre¬ 
ments  of  growth.  The  ten  values  of  Gj  for  the  smallest  a  are  represented 
by  a  horizontal  line  segment.  For  the  next  three  delamination  sizes,  both  G^ 
and  Gjj  were  constant  for  all  ten  increments  of  growth,  but  the 
ratios  Gj/Gjj  changed.  For  the  largest  delamination  size  modeled,  Gj 
and  Gjj  were  identical,  both  equal  to  one-half  the  total  G.  Hence, 
calculations  of  Gj  and  Gjj  varied  with  mesh  size. 
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EFFECT  OF  COUPLING  ON  STIFFNESS  LOSS 
The  stiffness  of  an  arbitrary  composite  laminate  may  be  calculated  from 
laminate  theory  [9,10]  as 


E 


1 

LAM  [A  '  ]t 
11 


(Bl) 


where 


[A']  =  [A]  1  +  [A]  1[B][D*]-1[B][A]-1 
[D*]  =  [D]  -  [B] [A]-1[B]  . 

and  [A],  [B],  [D]  are  the  extensional,  coupling,  and  bending  stiffness 

matrices,  respectively,  defined  in  reference  [9]«  If  the  laminate  is 
symmetric,  the  [b]  matrix  vanishes,  and  equation  (Bl)  reduces  to 
equation  (l). 

For  an  arbitrary  laminate  containing  delaminations  in  one  or  more 
interfaces ,  equation  (Bl)  may  be  used  in  equation  (2)  to  calculate  the 
stiffnesses  ,  ,  of  the  sublaminates  that  are  formed.  For  the 

[ ±30/±30/90/90 ]s  laminate,  with  delaminations  modeled  in  both  -30/90  inter¬ 
faces,  equation  (2)  becomes  equation  (3) 

8E(±30)2  +  3E(90), 

E*  =  - - - ^ 

11 


30 


APPENDIX  B 


If  the  :(±30)2  sublaminate  stiffnesses  are  calculated  using  equation  (Bl) 
and  material  properties  from  reference  [ll],  then  E*  =  0.69  E-y^,  which 
exceeds  the  value  of  E*  extrapolated  from  a  least-squares  regression  line 
for  the  stiffness  data.  However,  both  -30/90  interfaces  did  not  delaminate 
complete  (fig.  3).  As  shown  in  figure  20(a),  delaminations  shifted  from 
one  -30/90  interface  to  another  through  various  90°  ply  cracks.  Delamina,t ions 
were  believed  to  grow  in  this  manner  to  reduce  the  effect  of  bending-extension 
coupling  that  would  have  been  present  had  only  one  interface  been  cleanly 
delaminated  or  had  both  interfaces  delaminated  at  the  same  time.  Therefore, 
delaminations  were  modeled  over  both  -30/90  interfaces  (fig.  20(b))  but  the 
bending-extension  coupling  in  the  (±30)2  sublaminates  was  neglected. 

Hence,  the  stiffness  of  each  sublaminate  was  evaluated  using  equation  (l). 

This  resulted  in  E*  =  0.743  Eyyyj  which  agreed  with  the  value  of  E* 
extrapolated  from  a  least-squres  regression  line  for  the  stiffness  data 
(fig.  10). 

In  addition,  stiffness  loss  for  delaminations  of  various  size  was 
calculated  from  the  finite  element  analysis.  Concurrent  strip  delaminations 
were  assumed  in  both  -30/90  interfaces.  Uniform  axial  extension  was  also 
assumed  and,  hence,  bending-extension  coupling  was  ignored.  The  stiffness 
corresponding  to  a  particular  delamination  size  was  calculated  using  the 
technique  outlined  in  reference  [l6 ] .  Figure  21  shows  the  stiffness  loss 
predicted  by  the  linear  rule  of  mixutes  (R.O.M. )  equation  without  coupling 
(Ei  calculated  from  eq.  (l))  and  with  coupling  (Ei  calculated  from  eq.  (Bl)). 
The  finite  element  predictions  agreed  with  the  rule  of  mixtures  equation  that 
neglected  bending-extension  coupling.  Hence,  the  finite  element  results  also 
agreed  with  the  data  plotted  in  figure  10. 
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STIFFNESS  CHANGE  FROM  SECONDARY  MECHANISMS 
Ninety  Degree  Ply  Cracks 

The  spacings  of  90°  ply  cracks  in  the  laminated  portion  of  partially 
delaminated  [±30/±30/90/90]g  laminates  were  measured  directly  from  enlarged 
x-ray  photographs.  The  mean  spacing  is  plotted  in  figure  22  as  a  function 
of  delamination  size  (see  fig.  9)  for  three  quasi-static  test  specimens.  As 
the  delaminations  grew,  the  crack  spacings  decreased  to  within  values  pre¬ 
dicted  by  a  simple  one-dimensional  model  [6]. 

Near  the  90°  ply  cracks,  load  was  assumed  to  be  transferred  to  the 
neighboring  plies  and  then  back  into  the  90°  plies  to  form  the  next  crack. 

The  model  predicted  that  90°  ply  cracks  would  reach  a  characteristic  satura¬ 
tion  spacing,  K,  determined  by  the  laminate  stacking  sequence,  ply  thickness, 
and  ply  stiffness,  but  independent  of  load  history. 

For  the  [±30/±30/90/90]„  laminate,  the  model  predicted  a  saturation 

o 

spacing  of  0.838  mm  (0.033  in.).  For  a  [±30/130/903]  sublaminate,  repre¬ 
sentative  of  the  delaminated  portion  of  the  [±30/±30/90/90]  laminate,  the 
model  predicted  a  saturation  spacing  of  0.955  mm  (0.0367  in.). 

Stiffness  change  due  to  90°  ply  cracking  in  the  partially  delaminated 
specimen  was  estimated  hy  reducing  the  maximum  axial  stress  in  the  90°  plies 
based  on  the  density  of  cracks  present.  As  shown  in  figure  23,  the  transverse 
stiffness  of  the  90°  plies,  E22>  was  the  percentage  difference 

in  the  maximum  axial  stress,  calculated  from  the  one  dimensional  model,  at 
the  predicted  crack  spacing  K  and  at  one  half  the  spacing  K/2.  This 
difference  represents  the  reduction  in  load  carried  hy  the  90°  ply  once  the 
saturation  crack  pattern  has  formed. 
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For  the  [±30/±30/90/90]s  laminate,  a  twelve  percent  reduction  in  the 
axial  stress  in  the  90°  plies  was  predicted.  Reducing,  the  transverse 
stiffness,  E22,  of  the  90°  plies  by  twelve  percent  resulted  in  a  0.7  percent 
reduction  in  [±30/±30/90/90]s  laminate  stiffness.  Similarly  for  the 
[±30/±30/902]  sublsyninate,  a  l4  percent  reduction  in  the  axial  stress  in 
the  90°  plies  was  predicted.  Reducing  the  transverse  stiffness,  E22, 
of  the  90°  plies  by  l4  percent  resulted  in  a  2.2  percent  reduction  in 
the  [±30/±30/902]  sublaminate  stiffness. 

Because  the  estimated  contributions  of  90°  ply  cracks  to  stiffness 
reduction  in  the  [±30/±30/90/90]g  delamination  specimen  were  so  small, 
they  were  neglected  in  calculating  E*.  However,  for  other  laminates  and 
materials,  ply  cracking  may  have  a  greater  effect  on  axial  stiffness  and 
other  stiffness  parameters  [17,18]. 

+30/-30  delaminations 

The  largest  +30/-30  interface  delaminations,  a,  recorded  during  the 
four  quasi-static  tension  tests  were  measured  and  tabulated  in  table  2. 

Only  specimen  E20  had  a  +30/-30  delamination  of  significant  size.  Also 
shown  in  table  2  are  predicted  stiffness  values  for  the  [ ±30/1 30/90/90*] 

s 

laminate,  containing  both  -30/90  and  +30/-30  interface  delaminations.  The 
stiffness  was  predicted  using  equations  (2)  and  (4)  first  for  the 
+30/-30  delamination,  and  then  for  the  -30/90  delamination.  The  stiffness 
predicted  for  specimen  E20  was  three  percent  less  than  the  stiffness  pre¬ 
dicted  assuming  only  -30/90  delamination.  However,  the  data  agreed  best 
with  the  prediction  that  omitted  the  +30/-30  delamination.  Hence,  the 
effect  of  +30/-30  interface  delaminations  on  measured  [±30/130/90/90"] 

s 

stiffness  loss  was  neglected  in  calculating  E*. 
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TABLE  1.-  FATIGUE  GROWTH  RATE  MEASUREMENTS 
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TABLE  2.-  EFFECT  OF  +30/-30  DELAMINATION  ON  STIFFNESS  PREDICTION 
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Figure  5,-  Through-thickness  (7Z  distribution  at  edge  of 

[±30/ ±30/90/90] s  laminate. 


Figure  7,-  Rule  of  mixtures  analysis  of  stiffness  loss. 
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Delamination  size  and  stiffness  data  accamulatlon  during  quasi-static  loading 
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Stiffness  as  a  function  of  delaminat 
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function  of  delamination  size. 
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Figure  14,-  Typical  plot  of  delamination 
function  of  load  cycles. 


MEAN  OF  FOUR 


ro  2 
-o  -o 


5^ 


t. 


Figure  17,-  Effect  of  delamination  measurement  method  on  power  law  curve  fit, 
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Figure  18,-  Finite  element  discretization  of  specimen  cross  section. 
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Figure  21.-  Effect  of  coupling  on  stiffness  loss  prediction. 


[ ±30/ +30/ 90/ 90] .  LAMINATE 


:gree  ply  crack  spacing 
delamination  size. 


[±30/ ±30/90/90]  LAMINATE 


Figure  23.-  Estimate  of  stiffness  loss  from  90°  ply  cracks. 
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